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ABSTRACT 



> 

J^ ; We have observed SN 1998S which exploded in NGC 3877, with the Utrecht Echelle 

<^ • Spectrograph (6 — 7 km s^"^ FWHM) at the William Herschel Telescope and with 

the E230M echelle of the Space Telescope Imaging Spectrograph (8 km s^^ FWHM) 
Q\ . aboard the Hubble Space Telescope. Both data sets were obtained at two epochs, 

(~| , separated by 19 (optical) and seven days (UV data). 

*^| From our own Galaxy we detect interstellar absorption lines of Ca II K, Fe II 

O ■ AA2600, 2586, 2374, 2344, Mg I A2852, and probably Mn II A2576, at ulsr = -95 km s^^ 

'^ ' arising from the outer edge of the High Velocity Cloud Complex M. We derive gas-phase 

^ • abundances of [Fe/H] = —1.4 and [Mn/H] = —1.0, values which are very similar to 

warm disk clouds found in the local ISM. This is the first detection of manganese from 
^ . a Galactic HVC, and we believe that the derived gas-phase abundances argue against 

the HVC material having an extragalactic origin. 

At the velocity of NGC 3877 we detect interstehar Mg I A2852, Mn II 
AA2576, 2594, 2606, Ca II K and Na I D2,Di absorption lines, spanning a velocity 
range of —102 to +9 km s"^ from the systemic velocity of the galaxy (910 km s~^). 
Surprisingly, the component at —102 km s~^ is seen to increase by a factor of ^ 1 dex 
in iV(Na I) between 20-March-1998 and 8-April-1998, and in A^(Mg I) between 
4-April-1998 and ll-April-1998. 

Unusually, our data also show narrow Balmer, He I, and metastable UV Fe II 
P-Cygni profiles, with a narrow absorption component superimposed on the bottom of 
the profile's absorption trough. Both the broad and narrow components of the optical 



^Based in part on observations obtained with the NASA/ESA Hubble Space Telescope, obtained at STScI, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National 
Aeronautics and Space Administration, NAS5-26555. 



lines are seen to increase substantially in strength between March-20 and April-8. The 
broad absorption covers ~ 350 km s~^ and is seen in Mg II A2796, 2803 absorption 
as well, although there is no evidence of narrow Mg II emission forming a P-Cygni 
profile. There is some suggestion that this shelf has decreased in strength over seven 
days between April-4 and April-11. 

Most of the low-ionization absorption can be understood in terms of gas co-rotating 
with the disk of NGC 3877, providing the supernova is at the back of the disk as we 
observe it, and the H I disk is of a similar thickness to our own Galaxy. However, 
the variable component seen in all the other lines, and the accompanying emission 
which forms the classic P-Cygni profiles, most likely arise in slow moving circumstellar 
outflows originating from the red supergiant progenitor of SN 1998S. 

Subject headings: Galaxies: ISM — Galaxies: Individual (NGC 3877) — Stars: 
Supernovae: Individual (SN 1998S) — Stars: Circumstellar Matter — ISM: Clouds — 
Ultraviolet: ISM 



1. Introduction 

Supernovae which explode in galaxies other than our own Milky Way offer a unique 
opportunity to study in detail an extragalactic interstellar medium (ISM), since the supernova 
acts a bright background source against which absorption from gas along the sightline can be 
observed. Supernovae have been used successfully to probe galaxy ISMs both from the ground 
(e.g. [Wallerstein et al. 1972| ; [lenkins et al. 1984 P'Odorico et al. 1985| ; |Vidal-Madjar et al. 198^ ; 



[Steidel et al. 199C ; Meyer fc Roth 1991| ) using Ca II and Na I absorption lines, and more recently. 



from space with lUE ( Pettini et al. 1982| ; [Blades et al. 198q ). Since the launch of HST, we 



have been using the high resolution spectrographs to study the ISMs of nearby galaxies in more 
detail. Although a uniformity in the absorption line properties from an extragalactic ISM might 
be expected, considering that multiple sightlines through the Galactic ISM show very similar 
Mg II profiles, the characteristics observed have, in fact, been quite diverse. Towards SN 1992A, 
which exploded in the elliptical galaxy NGC 1380 in the Fornax Cluster, we detected no Mg II 
absorption to sensitive equivalent width limits ( [Bowen et al. 19"95| ), suggesting that the galaxy 
was devoid of neutral hydrogen, and that elliptical galaxies are unlikely to contribute significantly 
to the numbers of Mg II absorption systems seen at high redshift. Conversely, toward SN 1993J 
which arose in M81, we detected individual absorption components spanning a velocity range of 
~ 400 km s~^, grouped into 3 subcomplexes arising in gas from the disk of M81, from our own 



Galaxy, and from intergalactic tidal debris from interactions between M81 and M82 (Bowen et al 



1994). In this paper we add another sightline to our sample, that of SN 1998S, and once again 



find unexpected absorption line properties in the supernova's spectrum. 

The type II-L supernova SN 1998S exploded in the outer regions of NGC 3877, an edge-on Sc 



galaxy in the Ursa Major cluster with a heliocentric velocity of 910 it 6 km s ^ ( [Broeils fc van 



Woerden 1994). The SN sightline lies at a projected distance of 48 arcsec from the center of the 
galaxy, or 3.6 kpc, assuming a distance of 15.5 Mpc to NGC 3877 ( [Sanders fc Verheijen 1998 ). 



An image of the supernova and its host galaxy is shown in Figure |^. In this paper we present 
high-resolution UV and optical spectra of the supernova, taken to probe the ISM of the host 
galaxy. 




Fig. 1.— A 2 min R-band CCD image of NGC 3877 taken on 16-March-1998 by S. Bennett & A. 
Aragon-Salamanca at the Isaac Newton Telescope on La Palma. North-East is top-left. SN 1998S, 
lying along the spiral arm towards the South- West, is 48 arcsec from the center of the galaxy. 



2. Observations and Results 

Details of the observations are given in Table 1, which lists dates of observations, wavelength 
ranges covered, exposure times, and spectral resolutions. HST observations of SN 1998S were 
made at two epochs, each with the Space Telescope Imaging Spectrograph (STIS) and E230M 
echelle through the 0.2x0.06 slit using the NUV-MAMA detector. We used the orders extracted by 
the STSDAS pipeline calibration; however, the subtraction of the extracted background spectrum 
resulted in a poor signal-to-noise for the final spectra, since the background count rate was a 
significant fraction of the total count rate. For each order, we therefore replaced the background 



spectrum with a single value obtained from averaging the data values in the background spectra, 
and subtracted this from the gross extracted counts, thereby improving the signal-to-noise of the 
data. Where necessary, overlapping orders were summed after weighting pixels by the inverse of 
the noise derived from the pipeline extractions. Optical spectra were obtained with the Utrecht 
Echelle Spectrograph at the William Herschel Telescope: data acquired 20-March-1998 were taken 
with the E31 (31.6 lines mm~^) echelle and a 2.0 arcsec slit, centered at 5430 A; spectra taken 
8-April-1998 were made with the E79 grating (79.0 lines mm~^) and a 1.0 arcsec slit, centered 
at 5900 A and 4280 A. For both data sets, a SITel 2048x2048 CCD detector was used, and both 
were reduced with the NOAO echelle IRAF package.^ 

Tabic 1. Chronology of Observations 









AA 


Exp time 


FHWM 


T+ 


UT Date 


Obsv. 


Instrument 


(A) 


(min) 


(km s-i) 


(days)" 


20-Mar-1998 


WHT 


UES 


4190-6610 


20 


6.0 





4-Apr-1998 


HST 


STIS+E230M 


2277-3115 


195 


8.3 


+15 


8-Apr-1998 


WHT 


UES 


3620-5860 
4530-9050 


60 
60 


7.4 
7.4 


+19 


ll-Apr-1998 


HST 


STIS+E230M 


2277-3115 


349 


8.3 


+22 



" Days after maximum light, assumed to be 20-March-1998. Discovery was 3-March-1998. 



2.1. Cold interstellar gas in the Milky Way &; from HVC Complex M 

Fortuitously, at Galactic co-ordinates of I = 150.7, b = 66.0, the sightline to SN 1998S 



intercepts the group of High Velocity Clouds (HVCs) designated as MI in Complex M ([Wakker & 



van Woerden 199l| ). The complex is thought to lie at a height of < 3.5 kpc above the Galactic 



plane ( |Ryans et al. 1997 ; Danly et al. 1995 ; Keenan et al. 1995 ) and is the highest latitude 



conglomerate of HVCs known. In fact, the sightline passes through the very edge of the complex 
(as measured at 21 cm), but data from the Leiden/Dwingeloo 21 cm Survey ( [Hartmann fc Burton 



|1997 ) show that there exists strong emission at vlsr — ~87 km s ^, at least over the 36 arcmin 



beam centered at I = 150.5, b = 66.0. 

In our STIS spectra, we detect Fe II AA2600, 2586, 2374, 2344, Mg I A2852, and, probably, 
Mn II A2576 absorption lines at vlsr = —95 km s~^ (Figure ^). This velocity agrees well with 
the 21 cm emission from Complex M, which we reproduce at the bottom of Figure ^a. This is 
the first detection of manganese in an HVC ( Wakker & van Woerden 1997 ) . The Fe II A2382 line 



^The Image Reduction and Analysis Facility (IRAF) software is provided by the National Optical Astronomy 
Observatories (NOAO), which is operated by the Association of Universities for Research in Astronomy, Inc., under 
contract to the National Science Foundation. The Space Telescope Science Data Analysis System (STSDAS) is 
distributed by the Space Telescope Science Institute. 



from the Milky Way is also covered by our STIS spectra, but the HVC absorption is blended 
with redshifted Fe II A2374 from NGC 3877 and is disregarded in our analysis. Similarly, Mn II 
A2594 and Mn II A2606 absorption which might arise from the HVC are blended with other 
redshifted Fe II lines. Figure ^b shows that in the optical, weak Ca II K is also detected (at a 
statistically more significant level that the Mg I, since the signal-to-noise of the optical data is so 
much higher than the STIS data), but the Ca II H line is absent. Galactic Na I Di,D2 is badly 
contaminated by He I A5877 absorption and emission at the redshift of NGC 3877 (see §2.3), 
making it impossible to derive reliable measurements of absorption from the HVC. Equivalent 
widths, Wx, and errors, a{Wx), to the HVC lines are given in Table 2, derived from standard 
methods (e.g. Bowen et al. 1995 and refs. therein) — the Mn II A2576 is detected at exactly the 
3(T level. The Fe II lines from the HVC are heavily blended with lower velocity absorption, so their 
measured equivalent widths are given as lower limits in Table 2. 

As can be seen in Figure |2|, we also detect Fe II and possibly Mg I at i^lsr = +48 km s~^. 
These lines are highly unusual given that there is no 21 cm emission seen at this velocity. The 
absorption must arise either in gas with unusual ionization conditions or else must represent dense 
gas whose size is very much smaller than the width of the 21 cm beam, and whose contribution to 
the emission is therefore too weak to be detected. 

To estimate column densities of these clouds, we used our absorption line fitting routine 
MADRIGAL QBowen et al. 1995| ) to best fit theoretical Voigt profiles to Fe II AA2600, 2586, 2374, 2344, 



Mn II A2576, Mg I A2852, and Ca II H, K lines simultaneously. We used the sum of the April-4 
and April- 11 STIS data to improve the signal-to-noise of the UV lines, and adopted a Gaussian 
Line Spread Function (LSF) with a Doppler width of 6instr = 5.0 km s^^ for the E230M echelle, as 
described in the STIS Instrument Handbook. We assumed that all ions would have a single value 
of "J^LSR cind Doppler width b; the values of b for individual ions would be different if the widths 
of the lines are dominated by kinetic temperatures rather than bulk turbulent motions (since 
b = [{2kT/m) + 2f^]^'^, where m is the atomic mass of the species, T the kinetic temperature of 
the gas, and vt represents the turbulence in the gas). However, the lines are clearly resolved in 
our data, show no obvious substructure, and have b values which would not be expected for the 
low-ionization lines seen if representative of kinetic temperatures. It therefore seems likely that 
the width of the lines are dominated by turbulent motions, and we adopted a single b value for all 
species. 

The data are of insufficient quality to properly derive column densities for the absorption from 
the Milky Way near zero velocity; nevertheless, since absorption at these velocities contributes 
to the higher velocity lines, we need to include lower velocity components, even though the fits 
will not be precise. The results are given in Table 2, and fits to the Fe II, Mg I and Ca II lines 
are shown in Figure ^. The values of vlsr are indicated by tick marks. We find a velocity of 
■'^LSR = ~95 km s~^ for the HVC, and derive b = 12 A km s~^, which agrees well with a simple 
Gaussian fit to the 21 cm emission feature shown in Figure 0a, 621cm = 15-0 km s~^. For the 
cloud at positive velocities we measure fLSR = +48 km s~^, b = 10.0 km s~^. We have derived 
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Fig. 2. — Interstellar Fe II (a, left), Mg I, Mn II (b, right), and Ca II K, H absorption from the 
Milky Way along the sightline to SN 1998S. The values of ijlsr derived from Voigt profile fits are 
indicated by tick marks above the best fits themselves (solid line). The line labelled 'A' is Mn II 
A2594 arising in NGC 3877, and that labelled 'B' is Mn II A2576 also in the host galaxy. Absorption 
can be seen at a velocity of Vq = —101 km s~^ (f^LSR = ~95 km s~^), which arises from the 
edge of the High Velocity Cloud Complex M. The 21 cm emission along the sightline, centered at 
I = 150.5, b = 66, from the Leiden/Dwingeloo H I Survey (Hartmann & Burton 1997) is shown 
bottom of (a). 



errors in the column densities by re-fitting 500 synthetic hnes generated from the initial fit to the 
data, each spectrum given the same signal-to-noise as the original data (see Bowen et al. 1995). 
The errors quoted in Table 2 refer to the value of a derived from a Gaussian fit to the distribution 
of h and N measured from the synthetic data, and hence refer only to errors associated with the 
noise in the spectra. Errors in both h and t^lsr arising from these Poisson statistics are small, 
0.3 km s~^. In reality, the error in wlsr will be dominated by the wavelength calibration of the 
STIS spectrum, while the value of h will be affected by uncertainties in the STIS LSF. In the first 
case, the error in the dispersion solution and zero-point for the MAMAs is no more that ~ 1 pixel, 
or ~ 5 km s~^, while for the optical data, the residuals of the fits to the echelle arc lines suggests 
an error of only ~ 0.2 km s~^. In the second case, the accuracies of the STIS LSFs still remain 
largely unexplored. 

Table 2. Column densities towards Galactic Clouds 
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'^Measured values are lower limits since the lines are blended with components nearer zero km s ^ . 
Limits are 'i(j{Wx) 



With no detection of H I towards the cloud at wlsr = +48 km s~^, we can derive little 
about the physical conditions of the absorbing gas. However, for the absorption by the HVC in 
Complex M, we can make crude estimates of the metallicity of the HVC. The H I column density 
estimated from the Leiden/Dwingeloo 21 cm Survey is ~ 6 x 10^^ cm^^. Adopting the cosmic 
meteoritic abundance ratios log(Fe/H)0 = 



-4.49 and log(Mn/H)0 = -6.47 ( [Anders fc Grevesse 
1989), we derive gas-phase depletion values of Fe and Mn —1.44 and —1.03 respectively.^ In 
the Galactic interstellar medium both Mn and Fe are heavily depleted onto dust grains, and in 
fact, the values of [Fe/H] and [Mn/H] are remarkably similar to those seen in the local 'warm 
disk' clouds discussed by Savage & Sembach (1996). However, dust depletion is not the only 
explanation for sub-solar gas-phase depletions — a combination of low metallicity and less dust 
depletion within the HVC cloud could give similar values. Stellar abundance measurements toward 



'[X/H]=log(X/H)-log(X/H)c 



moderately metal-poor ([Fe/H] = —1.0) halo stars find Mn underabundant with respect to Fe 
with an average [Mn/Fe] = —0.48 (Ryan, Norris, & Beers 1996). High redshift measurements in 
damped Lya QSO absorption systems find [Mn/Fe] values consistently below zero, which has been 
interpreted as resulting from a combination of a warm Galactic disk-like dust depletion pattern 
superposed on a slightly sub-solar Mn/Fe abundance ratio ( Vladilo 199§| ). We present the ISM, 
stellar, damped Lya QSO [Mn/Fe] measurements, and our new HVC value, in Figure |3[ The 
derived [Fe/H] and [Mn/H] ratios give a gas-phase depletion value of [Mn/Fe] = 0.41 ±0.11 for the 
HVC, quite different from the values in metal-poor stars or extragalactic DLya systems. 
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Fig. 3. — [Mn/Fe] values in the Galactic ISM and damped Lya absorption systems plotted as a 
function of redshift. Our marginal Mn II detection toward the HVC in the SN 1998S sightline infers 
the [Mn/Fe] value indicated by a star symbol. The solar abundance ratio (Mn/Fe)© marked by 
a dashed line is from the meteoritic abundances of Anders & Grevesse (1989). Galactic halo star 
measurements are also shown, taken from Ryan, Norris, &: Beers (1996)| . Measurements of [Mn/Fe] 
for damped Lya QSO absorbers can be found in [Lu, Sargent, &: Barlow (1996), Lu et al. (1995)|, 



Meyer et al. (1995)| , Pettini et al. (1999)| , |Roth fc Songaila (1999)| , and |Steidel et al. (1995) . 



Although not conclusive, (since the point may only be an upper limit if the detection of Mn II 
is erroneous) our data are consistent with the HVC material at ^lsr = —100 km s^^ having a solar 
Mn/Fe abundance ratio and a dust depletion pattern similar to that observed in warm disk clouds 
in the Galactic ISM. This interpretation argues against the HVC material having an extragalactic 
origin as has been proposed for some HVCs (e.g. Blitz et al. 1999|). 



2.2. Cold interstellar gas in NGC3877 

The optical, and UV data in particular, show a wealth of absorption lines arising from 
NGC 3877, and in this section we focus on the low-ionization species detected. Figure ^ shows 
the optical Na I and Ca II K lines, with the abscissa showing the velocity of the components 
relative to the systemic velocity of the galaxy (taken to be 910 km s~^). Ca II H is clearly detected 
at the velocity of NGC 3877, but is contaminated by He emission &: absorption (see §2.3). The 
absorption spans a range of velocities covering a total spread of ~ 60 km s~^ for Na I and 
~ 100 km s~^ for Ca II, each complex clearly made up of a number of components. The range 
of absorption is consistent with the observed H I emission along the sightline, as deduced from 
the 21 cm emission maps of Broeils & van Woerden (1994), and marked by a bar in the bottom 
panel. We have plotted data taken at the two different epochs (Ca II was only observed April-8); 
the most striking thing about the difference between the two profiles is the increase in the most 
negative component of the Na I absorption over an interval of 19 days. This increase in absorption 
corresponds to an increase in a narrow absorption line at the bottom of the absorption trough of 
the He I A5875 P-Cygni profile. Although these profiles are discussed in more detail in §2.3, we 
introduce this feature, shown in the top panel of Figure ^, to highlight the unusual nature of the 
variable Na I component. 

In Figure ^3 we show absorption line profiles of Mg I A2852 and the Mn II triplet centered on 
Mn II A2594. We show only April-4's data for the sake of clarity, but there is evidence that the 
strength of the blue wing of the Mg I absorption has increased, which would correspond to the 
same increase seen in the Na I component. The same may true for the Mn II lines, but the increase 
in each of the multiplet lines is small, and within the variations from the low signal-to-noise of 
the data. We note too that Fe II AA2600, 2586 and Fe II AA2344, 2374, 2382 are also detected at 
the velocity of NGC 3877 in the STIS data. However, the lines are heavily blended with either 
zero-redshift Milky Way absorption from various species, or metastable Fe II absorption from the 
host galaxy. For these reasons we choose not to include the lines in our analysis. 



2.2.1. Derivation of column densities and Doppler parameters 

To derive column densities, N, and Doppler parameters, b, of the low-ionization gas, we 
apphed MADRIGAL to the Na I D2,Di and Ca II K hues taken 8-Apr-1998. A minimum of nine 
components were required. We initially attempted to fit both optical and UV absorption lines 
simultaneously. However, it was clear that the lower resolution and lower signal-to-noise of the 
STIS data meant that the spectra offered few data points with which to constrain the fits. Most 
of the information on the velocity structure in the extragalactic absorption is only evident in the 
optical data. We therefore first used a simultaneous fit to the Ca II and Na I complexes to derive 
velocities for individual components. Since the atomic masses of both calcium and sodium are 
different, in principle, the values of b derived from a fit need not be the same for each ion if the 
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Fig. 4. — a) (left) Plot of He I A5875, and interstellar Na I and Ca II K lines vs. v — tigystemic arising 
in NGC 3877 obtained with the UES at the WHT, taken between 20-Mar-1998 (heavy line) and 
8-Apr-1998 (light line). The Na I A5889 night sky line is marked by a ©. The position of Na I 
absorption from the Milky Way is also indicated in the top panel. A narrow Na I component at 
—102 km s^^ varies over 19 days; this is also seen in He I absorption, which varies over the same 
timescale. The development of a broad wing in He I absorption at v — ^systemic < —100 km s~^ is 
also apparent. Only one set of data covering Ca II were available, so variability in that species is 
undetected. The total extent of H I observed at 21 cm is represented by a bar in the bottom panel, 
while the dotted vertical line represents the systemic velocity of the host galaxy, b) (right) Plot of 
interstellar Mg I and Mn II arising in NGC 3877, taken from April-4's STIS data. 
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line widths are dominated by the kinetic temperature of the gas. Hence with vq constrained, we 
refit both the Ca II and Na I aUowing b k N to vary for most of the components. (Obviously, we 
assume that Vq is the same for each ion, although this assumption depends on the geometry and 
physical conditions in the cloud.) In fact, the Ca II components are poorly determined without 
the availability of the Ca II H line; a fit to both Ca II and Na I simultaneously, with b set to be 
the same for both species (but varying during the fit), merely results in nearly the same b values 
derived for Na I, along with adjusted column densities for Ca II. We thus re-considered the values 
of A^(Ca II) derived by using only the b values from Na I, but found little change in the resulting 
Ca II column densities. 

With plausible values of b and vq established for Na I and Ca II, we then fit the Mg I A2852 
and Mn II AA2576,2594 lines taken with STIS on Apr-11. Magnesium has an atomic mass next 
to sodium, so we kept Vq and b fixed to the values derived from the Na I fits — again, only the 
depth of the Mg I line is constrained in the line profile, and fixed b values are necessary to produce 
a meaningful fit. Manganese has an atomic mass closer to calcium, so we fit the Mn II line using 
the b values derived for the Ca II lines. The fits to the manganese lines have the advantage that 
we used two lines of the multiplet, which help constrain b and N, but the low resolution and poor 
signal-to-noise still make the fits uncertain. The fits are shown in Figure |5|a. 

N and b were then derived for Na I D2,Di data taken March-20. We found that we were 
able to well fit the data using exactly the same values of Vq, b, and N(Na I) derived from April's 
data, except for component 1, where A^(Na I) appears to have varied by nearly 1 dex in the space 
of 19 days. The profiles are shown in Figure gb. Finally, we refit 4-Aprirs Mg I again using the 
values derived from April's data but allowing N(M.g I) to vary for component 1. Formally, our fits 
suggest that A^(Mg I) has increased by a factor of ~ 200. In these data, however, the component 
is close to being saturated, and with the low signal-to-noise of the data, large increases in A^ can 
occur with very small variations in the noise of the line. Nevertheless, some constraint on the 
fit comes from the fact that the line has widened slightly as well as increased in depth. Hence 
although A^(Mg I) may be over-estimated, an increase over 1 — 2 dex still seems probable. All the 
derived physical parameters are given in Table 3. 

Ignoring the variable component (which contributes ~ 0.1 A), the equivalent widths of the 
optical lines seen in NGC 3877 are WA(Ca IIA3933)= 0.524±0.007, WA(Na IA5889)= 0.499±0.005, 
and WxCNa IA5895) = 0.319 ± 0.005 (Ca II H is blended with He from the host galaxy). The Na I 
equivalent widths are entirely consistent with the distribution measured along sightlines to O- and 
B-type stars in the Galaxy ( Sembach et al. 1993| ) while the Ca II equivalent widths are slightly 



higher than the median of the equivalent width distribution seen along extragalactic sightlines 
through the Milky Way [0.1—0.2 A for Ca II ( Bowen 1991| )]. Similarly, the total column densities 



measured in NGC 3877 [SlogiV(Na I) = 12.6, Slog A^(Ca II) = 12.8], again excluding the variable 
component, are close to those measured through the entire length of the Galactic plane ( pembach 



& Danks 1994). In this sense, the optical lines suggest that the ISM of NGC 3877 along the line 



of sight is similar to our own. 
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Fig. 5. — a) (left) Fits of theoretical Voigt profiles to S-April's optical and ll-April's UV data. 
The solid line shows the resulting line profile from a convolution of 9 components. The lighter 
lines show the profile of each component for the values given in Table 3, to highlight how each 
contributes to the final convolved profile, b) (right) Fits of theoretical Voigt profiles to 20-March's 
optical data. The values of ^Na i and iV(Na I) are the same as 8-April's data shown in (a) except 
for component 1, which has clearly increased in column between the two epochs. The dotted line 
represents S-April's fit to component 1 to highlight the change. 



Table 3. Column densities from absorption lines in NGC 3877 
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logN 
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b 
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b 
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logAf 


No. 


(km s-i) 


(km s-i) 


(km s-i) 


logN 


3<7(N) 


(km s-i) 


(km s-i) 




logN 


3<7(N) 


Mgl 


3<7{N) 


Mnll 


3<7{N) 


Mgl 


Nal 


1 


-102.1 


2.7 


0.2 


11.89 


0.04 


3.7 


0.2 




12.00 


0.06 


13.80 


0.39 


13.20 


1.56 


11.5 


11.0 


all epochs 


2 


-90.1 


7.6 


1.5 


11.64 


0.17 


b 






11.96 


0.19 


12.36 


0.57 


12.96 


0.21 






3 


-78.7 


15.7 


0.7 


12.09 


0.06 


b 






11.99 


0.28 


12.28 


1.37 


13.03 


0.08 






4 


-71.3 


1.3 


0.7 


11.90 


0.09 


3.0 


0.6 




12.06 
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0.40 
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0.28 


<12.4 






CO 
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11.88 
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0.21 


12.01 


0.23 
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11.90 
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12.07 


0.14 


<12.0 
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0.4 




11.79 


0.10 


11.90 


0.15 


< 12.0 









^Velocity relative to systemic, taken to be 910 km s ^ 
^6(Na I) used to constrain fit 
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In Figure | we show the A^(Mg l)/iV(Na I), A^(Ca II)/7V(Na I), and N{Ca II)/7V(Mg I) ratios 
from the fits in Table 3. The use of the first of these, A^(Mg l)/A^(Na I), as a diagnostic of gas 
temperature has been discussed by Pettini et al. (1977)| , who showed that low ratios in the Galactic 



disk could be explained if absorption occurred only in cool (T < 80 K) H I clouds. The ratios 
found towards SN 1998S show a mixture of values, suggesting both cool quiescent clouds and 
warmer halo structures. Similarly, the A^(Ca Il)/A^(Na I) ratio, usually interpreted as an indicator 
of the calcium gas-phase abundance, also shows a wide range; both, however, suggest a difference 
in physical conditions above and below ~ —50 km s~^. We consider these ratios more fully below. 



2.2.2. Interpretation 
At the position of the SN sightline, the observed H I from 21 cm emission in NGC 3877 is 



measured from 730 — 910 km s ^ ( Broeils fc van Woerden 1994 ), or Av = v — ^systemic of —180 



to km s~^. The galaxy is highly inclined, i = 84°, yet the distance of the SN from the major 
axis is still non-negligible, which suggests that the true radial distance, R, of the object from the 
center of the galaxy is large, as much as 10.3 kpc, assuming it lies in the galactic plane. If so, the 



SN is near the edge of the H I radius of the galaxy, measured to be ~ 10 kpc (Sanders &: Verheijen 



1998). The question remains, however, whether the SN sits at the front or back of the disk, and 
whether the low-ionization absorption can be understood in terms of disk-gas along the line of 
sight. If the true distance of the SN in the disk is 10 kpc, then the de-projected angle between it 
and the major axis of the galaxy is ~ 71°. The true radial velocity, v{R), of gas in the disk of 



NGC 3877 at 10 kpc is 175 km s ^, as measured from the rotation curve ( Sanders fc Verheijen 



1998]) . Hence the projected velocity, Vq, along the line of sight of gas close to the supernova's 
environment would be Vq — iJays = v{R) sin i cos ~ —56 km s~^, where 9 is the angle between a 
point in the disk along the sightline and the center of the disk. This is true whether the SN is in 
the edge of the disk closest or furthest from us, and the velocity lines up well with the red edge 
of the Na I complex in Figure |^. Perhaps the simplest argument to suggest that the SN must lie 
on the trailing side of the galaxy is to point out that at 10 kpc, the rotation curve is flat, and 
were the SN to lie on the leading edge, all the disk gas which lies between the probe and us would 
essentially be at the same velocity, ~ —56 km s~^. From Figure H, this does not appear to be 
the case. If the galaxy were exactly edge-on to us, and the SN occurred at the back of the disk, 
we might expect to see absorption from gas through the entire disk. The most negative velocity 
gas would occur where the sightline passes closest to the center of the galaxy (i.e. 9 = 0). This 
value (essentially the projected distance of the SN from galaxy center) is 3.5 kpc , and v{R) is 
still ~ 150 km s^^. The projected velocity of the gas would then be —146 km s~^ from systemic, 
and we do not see any absorption at these velocities. In fact, given that the galaxy is not edge-on, 
and assuming that the H I disk has a height similar to that in our own Galaxy, ~ 1 kpc, the 
sightline never intercepts gas near the center of the galaxy — it probes the ISM over perhaps only 
half a disk radius, depending on where exactly the supernova sits at the back of the disk. Then 
for R = 5 kpc, and v{R) = 160 km s~^, Vq — i^sys = —90 km s~^, which is, remarkably, the same 
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Fig. 6. — Plots of the derived column density ratios for the 9 low-ionization components seen in 
NGC 3877. The dotted lines in the top figure represent the range of values of A^(Mg I)/A^(Na I) seen 
by Pettini et al. (1977)| towards Galactic stars. The dashed lines in the middle figure represent the 
range of A^(Ca II)/A^(Na I) seen in the Milky Way for absorbing components at peculiar velocities 
of Av < 50 km s~^ ( Sembach &: Danks 1994 ). The corresponding component number given in 
Table 3 is shown top of upper panel. 
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as the most blueward non-variable component we see. 



If this simple model of absorption by the back quarter or so of the galactic disk is correct, 
then it also explains to a large degree the distribution of column density ratios vs. velocity as 
plotted in Figure ^. The components between —40 and —90 km s~^ (i.e. numbers 2 — 4 and 
6-7) all lie in the disk of the galaxy and have N{Ca U)/N(Na I) and N(Mg I)/A^(Na I) ratios 
similar to those seen in the disk of the Milky Way. Two other components, numbers 8 and 9, 
have velocities which do not fit in with the rotation of the galaxy — these components have high 
peculiar velocities, relative to the gas in the disk, similar to those seen in our own galaxy, and have 
correspondingly high A^(Ca Il)/A^(Na I) and A^(Mg l)/A^(Na I) values. Component 5 would then 
be interpreted perhaps as arising from a warmer intercloud medium. This just leaves us to explain 
the origin of the variable component 1, which has no apparent relation to the co-rotating ISM gas 
in NGC 3877. This component corresponds to many other unusual features, however, which we 
now discuss. 



2.3. Circumstellar absorption and emission from SN 1998S 

Although the initial aim of our observations was to examine the UV absorption from the 
interstellar medium of NGC 3877, our data have revealed unusual absorption line profiles which 
have more to do with the environment of the supernova itself. In Figure ^ we show the Mg II 
absorption both from our own Galaxy and from NGC 3877. The system is complicated by Mg II 
A2796 at the velocity of NGC 3877 being contaminated with Mg II A2803 from our Galaxy at rest. 
To highlight the possible degree of contamination, we have shifted and overplotted the Galactic 
Mg II A2796 profile at the position expected for Galactic Mg II A2803, as shown by the hatched 
region in the top panel. Hence the absorption seen in this region is a superposition of saturated 
lines from both Mg II A2796 absorption in NGC 3877 and Mg II A2803 from the Milky Way. This 
same blend of absorption is also seen in the middle panel, with Mg II A2803 from the Galaxy now 
at rest, but contaminated with Mg II A2796 from NGC 3877. 

Despite this confusion, it is clear that there exists a 'shelf of absorption spanning 
~ 350 km s~^ in addition to the saturated components expected from the ISM of NGC 3877 itself 
at Vq = 800 — 900 km s~^. We see a smooth decline in the strength of the absorption as the 
velocity from systemic increases. There is also some evidence that the strength of the absorption 
in this region has declined between April-4 and April- 11, although we note that defining the 
continuum in regions of wide absorption is difficult. 

Features in the optical data also show unusual conditions. At the bottom of Figure [^ we 
overplot the Ha profile taken April-8th. The 'shelf seen in Mg II AA2796, 2803 can also be seen in 
Ha absorption, but now part of a classic P-Cygni profile defined by Ha emission to the blue. This 
absorption can be seen in other Balmer lines too, in Ha, H/J, H7 and US. The absorption lines 
do not result from incorrect background subtraction of emission from other parts of the galaxy 
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Fig. 7. — Upper two panels show the Mg II hne profiles obtained with STIS on April-4 and April-11. 
See text for an explanation of the hatched areas. Bottom panel shows the optical spectra covering 
the emission and absorption from Ha. Mg II A2791.60 (labelled '1') and Mg II AA2798.754, 2798.823 
('2') are also seen, arising at the same velocity as the variable, narrow Na I component No. 1 shown 
in Figure |5[ 
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filling the slit — extraction of the supernova signal without background subtraction also shows the 
P-Cygni profiles. We show these profiles in Figure |8|, with the data from the two different epochs 
plotted over each other. (He is also present, but blended with Ca II H from the ISM of NGC 3877, 
so is not shown in Figure ^.) There is no doubt that the absorption complex has strengthened 
significantly between March-20 and April-8, in contrast to the Mg II which has declined. This 
suggests that the gas forming the Mg II has increased in ionization, with more of the magnesium 
forming Mg III. 

The exact level of Ha taken March-20 is uncertain, because no flux calibrators are available, 
and there is difficulty in defining where best to fit the continuum. There is, however, a possibility 
that the Ha shows weak excess emission as a wide red wing covering ~ 600 km s~^ on March-20, 
emission which has disappeared by April-8. These P-Cygni profiles, so well displayed in Balmer 
lines, can also be seen in many He I lines (e.g. He I A5875 in Figure ^); we list all the identified 
emission lines in April-8's data in Table 4, and note whether the lines show the P-Cygni profile. 
Some do not, and these emission lines are likely to be due to either the H II region in which the 
supernova exploded (see below), or from other emission line regions along the line of sight not 
directly connected with the supernova's environment. 

What is also clear from these P-Cygni profiles is that there exists strong, narrow absorption in 
addition to the weaker, wider shelf already described, located at the bottom of the P-Cygni trough. 
Although the exact center of this absorption is hard to measure, because it is superimposed on 
the emission line, it seems clear that it lies at the same velocity as the variable Na I component 
at Vo — Wsys = —102 km s~^ discussed in §2.2. This can best be seen in Figure ^ where the sharp 
absorption in He I A5875 is coincident with the varying Na I component. Further, there are other 
absorption lines in our HST spectra which coincide with this varying component. For example, 
there are the rarely seen Mg II A279I.60 and Mg II AA2798.754, 2798.823 lines (labelled '1' and '2' 
in Figure |^ which arise at this same velocity. 

Even more unusual are the many uvl, uv2, uv62, uv63, uvI58 and uvl59 multiplet Fe II lines 
seen in the STIS spectra, also with P-Cygni profiles. All the transitions we detect arise from the 
a'^H, a^D, and a'^D states. Examples of these lines are shown in Figure 0, where we have plotted 
all the lines of a given multiplet shifted from rest to a velocity of 850 km s~^. Both the broad 
shelf and the narrow absorption is seen in these profiles, as detected in the other species. There 
is no clear difference between the profiles between April-4 and April-Il. Although the data are of 
too low a signal-to-noise, and the continuum of the P-Cygni profiles too ill-defined, to measure 
physical parameters accurately, the narrow component is clearly resolved in the data. Crudely, the 
lines appear to be a factor of 2 wider than the instrumental resolution, which gives a lower limit 
of 6 > 6 km s~^, assuming that the spectra have a resolution of 1.7 pixels FWHM, as given by 
the STIS Instrument Handbook. This is much wider than that deduced from the Na I component 
alone, {b ~ 2.7 it 0.2 km s^^, Table 3) which suggests a much more complicated structure to the 
narrow component than a single absorbing cloud. 
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Fig. 8. — Sequence of Baliner lines observed 19 days apart with the UES at the WHT. 
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Fig. 9. — Two figures showing some of the metastable Fe II hnes detected in April-4's STIS spectra. 
All the lines arising from the uv62, uv63 (top) and uvl (bottom) transitions are plotted over the 
wavelength ranges shown, shifted by 850 km s~^, regardless of whether we consider them detected. 
Also shown are the multiplet number and the rest wavelength of the line. The uv62 and uv63 lines 
involve a change from a a'^D to a z^F° state, while the uvl lines arise from a change in energy 
from a^D to z^D". 
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Table 4. Emission lines identified 
in 8-Apr-1998 optical data 



Line 


P-Cygni 


[Ne III] A3867 


n 


He I A3870 


n 


He I A3887 + H 8 


y 


He I A3963 


y 


He 


y 


He I A4025 


y 


H5, H7 


y 


[0 III] A4361 


n 


He I A4470 


y 


He I A4711 


y 


He I A4290 


? 


[Fe HI] A4656 


n 


up 


y 


[0 III] A4957 


n 


[0 III] A5005 


n 


He I A5014 


y 


[Fe HI] A5268 


n 


[N II] A5753 


n 


He I A5889 


y 


Ha 


y 


[N II] A6581: 


n 


I A7769 


n 



It seems clear that all these features — the emission lines and the narrow, variable absorption 
component superimposed on a broader shelf of absorption — are unlikely to have anything to do 
with the supernova ejecta itself. Typically, ejecta material is moving at ~ 10^ km s~^ at the epochs 
at which we observed SN 1998S, hence the P-Cygni profiles would be considerably wider than the 
~ 30 km s^^ FWHM observed here. Indeed, the far-UV G140L spectra show the expected broad 
P-Cygni profiles from the expanding photosphere of the supernova in a variety of ions. 

Instead, it seems likely that the P-Cygni profiles arise from circumstellar material, blown 
off from the progenitor star. Such an interpretation is not unprecedented. Narrow Balmer lines 
and He I A5876 have been seen before in a supernova spectrum, towards SN 1984E in NGC 3169 
(Dopita et al. 1984). Their detection was cited as evidence of a superwind flowing out from the 



precursor star. In this case, the wind had a velocity of some 3000 km s~^, much greater than we 
observe towards SN 1998S. UV emission lines towards SN 1979C in M 100 were also interpreted 
as being due to a UV-emitting shell of compressed gas from pre-existing circumstellar material 



(Panagia et al. 1980), and indeed, a narrow Ha P-Cygni profile may also have been seen (Branch 



et al. 1981 ). Again, however, the width of these UV emission lines ranged from being unresolved 
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at a limit of 1000 km s to a maximum of 4000 km s , larger than we see towards SN 1998S. 

The circmxLstellar absorption detected towards SN 1998S is obviously far more quiescent than 
these superwinds. Stellar winds from, e.g., O- & B-type stars, or Population I Wolf-Rayet stars, are 
known to have high mass-loss rates and strong outflows, giving rise to the classic P-Cygni profiles. 
However, these outflows move with a speed of a few xlO^ km s~^ (see, e.g. Lozinskaya 1992), 



again, much higher than we observe towards SN 1998S. Perhaps the most interesting comparison is 



the detection of the same Fe II UV multiplets seen in the winds of A- type supergiants ( Talavera fc 
Gomez de Castro 1987 ; Verdugo et al. 1999| ) . Although not all A-type stars show these absorption 
lines, those that do are characterized by outflows moving at only a few hundred km s^^, along 
with variability in the absorption. It seems likely that this type of outflow is responsible for the 
complicated emission and absorption seen in NGC 3877. 

Another possibility is that at least some component of the absorption arises from a 
proto-planetary disk surrounding the progenitor. Observations towards /3-Pictoris and other 
A-type and late 0-type stars show the same Fe II metastable lines as we find towards SN 1998S, 
as well as variable absorption (e.g. Lagrange et al. 1996, Kondo &: Bruhweiler 198^ , Cheng et al, 



1997). If so, the absorption we see might well mark the presence of an early solar system shortly 
before its destruction by the supernova. The most obvious argument against this interpretation 
is simply that the absorbing shell or wind is probably much further away from the progenitor 
than would be expected from a proto-planetary disk. The existence of the Fe II lines in our last 
STIS spectra taken April-11 means that the ejecta from the supernova had not yet reached the 
absorbing material, some 39 days after the supernova was detected. Assuming simply that the 
ejecta moved with a constant velocity after outbreak, at, e.g., 5,000 km s~^, then for the absorbing 
material to remain undisturbed (i.e, assuming that the ejecta would have catastrophic effects on 
the absorbing material), the distance of the material from the progenitor would be > 100 AU. 
This is much larger than the few AU postulated for the /3-Pictoris disk. 

We also note that the shell has a velocity relative to the progenitor in our simple model of 
where the supernova sits in the galaxy, discussed in §2.2.2. We associated the progenitor to be at 
a radius in the galaxy where the ISM gas has a projected velocity along the line of sight to be 
V — ^systemic ~ — 50 km s~^] if this gas forms a reference frame for the progenitor, then the variable 
absorption occurring at f — ^systemic ~ —100 km s^^ is ~ 50 km s~^ different, a value which 
could therefore be taken to be the shell's speed. This, combined with the broader absorption 
shelf, would suggest a dense shell with a speed of about 50 km s^^ and a more highly ionized shell 
moving at ~ 300 km s~^. 

There are obvious complications in trying to interpret the observed data with stellar outflows. 
First, the material does not necessarily form a single bubble of gas; towards SN 1987A, for 
example, it is believed that the slowly expanding wind from the red supergiant phase of the 
progenitor is being caught by a faster wind from its blue supergiant phase, producing two shock 
fronts. There is also the obvious fact that the wind must hit the ambient ISM surrounding the 
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progenitor, producing another type of shock front. These process are not weh understood in our 
own Galaxy, and applying them to the features observed towards SN 1998S is clearly not easy. 
Second, towards a supernova, the ionization structure of the gas is complicated by the addition 
of the UV flash occurring at the initial explosion of the supernova. The effects of this ionizing 
source are superimposed on the ionization of the winds and shells caused by the progenitor and 
the earlier mechanical interactions. Third, in attempting to link the variability of the narrow 
absorption with the types of variations seen in the outflows of type-A supergiants, there is the 
added complication that the UV source against which the absorption is seen is expanding. This 
could in principle produce — apparently — variable absorption features, depending on the density 
structure within the outflow (whether clumps exist in the outflow) and the covering factors of 
the clumps. Finally, although the Fe II lines can be found in the spectra of A-type stars and 
their circumstellar material, they rarely show the well defined P-Cygni profiles we see towards 
SN 1998S. 

Hence, although we believe that we are seeing unique signatures of the outflows from the 
supergiant progenitor of SN 1998S, the exact details of the structure of the circumstellar material 
are beyond the scope of this paper. 



3. Summary 

We have observed SN 1998S which exploded in NGC 3877 with the UES (6-7 km s'^ 
FWHM) at the WHT and with the E230M echelle of the STIS (8 km s"^ FWHM) aboard HST. 
Both data sets were obtained at two epochs, separated by 19 (optical) and seven days {HST data) 
(see Table 1). We summarize our results as follows. From our own Galaxy we detect: 

• interstellar Ca II K, Fe II AA2600, 2586, 2374, 2344, Mg I A2852, and probably Mn II A2576, 
absorption lines at vlsr = —95 km s~^ arising from the outer edge of the High Velocity 
Cloud Complex M. We derive gas-phase abundances of [Fe/H] = —1.4 and [Mn/H] = —1.0, 
values which are very similar to warm disk clouds found in the local ISM. 

This is the first detection of manganese from a Galactic HVC, and we believe that the values 
of the iron and manganese gas-phase abundances argue against the HVC material having an 
extragalactic origin. 

At the velocity of NGC 3877 we detect: 

• interstellar Mg I A2852, Mn II AA2576, 2594, 2606, Ca II K and Na I D2,Di absorption lines, 
spanning a velocity range of —102 to +9 km s~^ from the systemic velocity of the galaxy 
(910 km s^^). In particular, the component at —102 km s^^ is seen to increase by a factor 
of ^ 1 dex in A^(Na I) between March-20 and April-8, and in A^(Mg I) between April-4 and 
April-11; 
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• Balmer and He I P-Cygni profiles in the optical, with a narrower absorption component 
superimposed at the bottom of the absorption trough of the profile. Both the broad and 
narrow components are seen to vary substantially in their absorption strength between 
March-20 and April-8; 

• a broad shelf of Mg II A2796, 2803 absorption spanning a velocity range of ~ 350 km s~^, 
with an apparent optical depth gently declining towards more negative velocities. This shelf 
covers the same velocity extent as the Balmer and He I lines in the optical (which form the 
blue side of a P-Cygni profiles) but there is no evidence of any narrow emission forming a 
Mg II P-Cygni profile. There is some suggestion that this shelf has decreased in strength over 
seven days between April-4 and April- 11; 

• metastable Fe II absorption lines from the UV multiplets uvl, uv2, uv62, uv63, uvl58 
and uvl59. These too show the same P-Cygni profiles seen in the optical, and are again 
superimposed with the narrow component at the trough of the broad absorption part of 
the profile. These lines do not obviously vary over the seven days between the two STIS 
observations; 

• rarely seen Mg II A2791 and Mg II AA2798.75, 2798.82 hues. These show P-Cygni profiles on 
April- 11 which do not exist on April-4. 

Most of the low-ionization absorption (the Ca II, Na I, Mg I, Mn II lines) can be understood in 
terms of gas co-rotating with the disk of NGC 3877, providing the supernova is at the back of the 
disk as we observe it, and the H I disk is a few kpc thick. However, the variable component seen in 
all the other lines, and the accompanying emission which goes to form the classic P-Cygni profiles, 
most likely arises in slow moving circumstellar outflows originating from the red supergiant 
progenitor of SN 1998S. 

We would like to thank Steve Smartt for taking the UES Service data on April-8, and in 
particular Sandi Catalan and Janet Wood for using their UES observing time on March-20. 
Thanks also to Peter Garnavich for sharing initial low resolution observations, Steve Bennett & 
Alfonso Ar agon- Salamanca for measuring SN magnitudes shortly after detection, Stefi Baum for 
sheparding the Phase 2 program at STScI, and the referee, Steven Federman, for a careful reading 
of our paper. The Fe II P-Cygni profiles listed in Table 4 were initially identified using Peter van 
Hoof's Atomic Line List v2. 02 web page (www.pa.uky.edu/ ""peter/atomic/index. html). 
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